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We characterize the phenomenon of “crowding” near the largest eigen¬ 
value Amax of random NxN matrices belonging to the Gaussian /3-ensemble 
of random matrix theory, including in particular the Gaussian orthogonal 
(/3 = 1), unitary (/3 = 2) and symplectic (/3 = 4) ensembles. We fo¬ 
cus on two distinct quantities: (i) the density of states (DOS) near Amax, 
PDOs(r, N), which is the average density of eigenvalues located at a distance 
r from Amax (or the density of eigenvalues seen from Amax) and (ii) the prob¬ 
ability density function of the gap between the first two largest eigenvalues, 
Pgap (r, N). Using heuristic arguments as well as well numerical simula¬ 
tions, we generalize our recent exact analytical study of the Hermitian case 
(corresponding to /3 = 2). We also discuss some applications of these two 
quantities to statistical physics models. 

PAGS numbers: 02.10.Yn, 05.40.-a 

1. Introduction 

During the last 20 years, there has been an important activity, both in 
mathematics and in physics, aiming at describing the fluctuations of the 
largest eigenvalue in ensembles of random matrices m El El nisi El E]. The 
most studied ones in this context, which we also focus on here, are probably 
the Gaussian /3-ensembles where the joint probability density function (PDF) 
of the N real eigenvalues Ai, • • • , Aat is given by: 

-Pjoint(Ai,A2,...,AAr) = ^P[|Ai - Ajl^ exp^-^y^A^^ , (1) 

where the normalization constant is Zfq = + 

fi/2)^ n^i + /^3/2) and where /3 > 0 is the Dyson index that can take 

* for Random Matrix Theory : Foundations and Applications, Krakow. 


( 1 ) 



2 


any real value. The classical values correspond to /3 = 1, 2 and 4, associated 
respectively to the Gaussian Orthogonal Ensemble (GOE), the Gaussian 
Unitary Ensemble (GUE) and the Gaussian Symplectic Ensemble (GSE). 
Note that for arbitrary /3, it is possible to associate a matrix model to Q 
(namely tridiagonal random matrices introduced in [8]). The fluctuations of 
the largest eigenvalues Amax = niaxi<j<jv A*, characterized by its cumulative 
distribution Fi\f{y) = Proba.[Amax < y], are now well understood. Indeed, 
we have now a precise characterization of the typical fluctuations of Amaxi 
when |Amax — V2A^| is of order which are described the TW 

distributions mm as well as of the large deviations of Amaxj when [Amax ~ 
V2N\ is of order 0{y/N), where F^iy) is described by large deviations 
functions [3l IH |5l [H [7] (both left and right tails). 



Figure 1. Different quantities characterizing the "crowding" near the largest eigen¬ 
value Amax studied in this paper: (i) the mean density of states /Odos(^iN) 
such that pDos{r, N)dr is the mean number of eigenvalues located in the inter¬ 
val [Amax — r — dr, Amax — r] and (ii) the PDF PGAp(r, N) of the spacing between 
the two largest eigenvalues, PGAp{r, N)dr = Pr .[(Ai^^r — A 2 , 7 v) G [r, r -|- dr]]. 


However in several situations, some observables related to such spectral 
statistics Q, might be sensitive not only to Amax but also to near-extreme 
eigenvalues, whose amplitudes are smaller but close to this largest eigenvalue. 
This general feature, not restricted to random matrices, has led physicists 
to study the phenomenon of “crowding” close to the maximum |9]. This was 
studied in detail for the case of independent and identical random variables 
|9] and more recently for Brownian motion |TU] . In Ref. [11] we proposed to 
characterize this phenomenon of “crowding” in random matrix models, by 
studying the density of near extreme eigenvalues, also called the density of 
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states, /3 dos(^)-^)) defined as (HUH] (see Fig. 

1 ^ 

pDOs{r, N) = ^ (-5(Amax “A* - r)) , (2) 

i=l 

where imax is such that = Amax and (• • •) means an average taken with 

the weight in Q. It is normalized according to 

dr pDosir, N) = 1 . (3) 

Incidentally, we showed that the DOS is related to another interesting 
observable characterizing also the crowding to Amaxi namely the PDF of 
the gap between the two largest eigenvalues m Let us denote by Amax = 
^i,N > ^ 2 ,N > • • • > ^N,N and by di^N = Ai^tv — ^ 2 ,n the first gap (see Fig. 
[^. Its PDF is denoted by pQApir, N), such that Pr.[di^jv £ [r, r + dr]] = 
PGApif, N)dr. It is then possible to show the following relation: 

PGAp{r, N) = {N - l)pDOs{-r, N) , (4) 

and we refer the reader to Ref. m for the derivation of this relation. 

In Ref. m, we focused essentially on the case of GUE (/3 = 2) where we 
computed exactly these two quantities /9 dos(^ 5-A) and pGAp(?’j-At) both for 
finite N and in the large N limit, using a method based on (semi-classical) 
orthogonal polynomials, which were previously introduced in Ref. m- Here, 
we generalize these results to the case of generic (3, which we mainly study, 
for large N, using heuristic arguments, that are confirmed by numerical 
simulations (for /3 = 1,2 and /3 = 4). Before presenting our results we 
discuss two applications of these two near extreme observables 

2. Two applications of Pdos(^? AI) and Pgap(^? AI) 

As we discuss it here, /Odos(^) -A) and PGApii^, 3V) play an important role 
in the minimization of quadratic forms on the sphere |13| and in the study 
of a spherical mean-field spin-glass m respectively. 

2.1. Minimizing a quadratic form on the N-dimensional sphere 

Let us consider the problem which consists in minimizing a random 
quadratic form on the A^-dimensional sphere Sn : 

^ N N 

iL[s] = — — JijSiSj , s = 'Y, ^ 1—^1 
i,j=l i=l 



( 5 ) 
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where J is a matrix belonging to the Gaussian /3-ensemble (with /3 = 1,2 
being the more natural values in this case) - and here we choose the matrix 
elements Jij of zero mean and variance of order 0{1/N) To take into 
account the spherical constraint (|^ we introduce a Lagrange multiplier z 
such that we have to minimize 

1 ^ 

H[s, z] = Ji,jSiSj z (s ^ - N) (6) 

with respect to sand z. It is then straightforward to show that 

inin H[s, z] = #[Smax, Zmax] = -N-^ , (7) 

s ,z Z 

where Smax and ^max are such that 


J ^max — '^max^max ; 

_ -^max 

-^^max — n • 


( 8 ) 


If we now look at the Hessian matrix 6^H/6siSsj evaluated at the minimum 
Smax) Zmaxi One Can show that its spectrum is given by 


Sp 


d^H 


SsjSs^ 


’max 5'^max 


— {0, Ajnax Ai, Aiuax A 2 , ' ' ' , Ajjiax Ajv} • (9) 


Hence the DOS, /9 dos(^)-A^) in Eh- @ is the average density of eigenvalues 
of the Hessian evaluated at the minimum (except the trivial eigenvalue 0). 
It is thus natural to expect that the DOS plays an important role in the 
relaxational dynamics of a system driven by such a quadratic form Q , which 
corresponds precisely to the fully connected p-spin spherical spin-glass model 
with p = 2 (where the Sj’s correspond to spin variables coupled to each other 
via the matrix J) [laiHini- 


2.2. Overlap distribution in the fully eonnected p = 2-spherical spin glass 

model 


In Ref. |14) . the authors studied the equilibrium properties of the p = 
2-spherical spin glass model (or spherical Sherrington-Kirkpatrick model) 
described by the partition function associated to the above Hamiltonian ([^ 



^ note that the joint PDF Q corresponds instead to matrix elements of order 0(1). 









where the 6 function ensures the spherical constraint ([^. They focused on 
the overlap Q, which is the order parameter characterizing the spin-glass 
order, defined by 


N 

2=1 


( 11 ) 


where and represent the spins at site i in two distinct equilibrium 
conhgurations with the same realization of the couplings Jij. In Ref. |l4| it 
was shown that, for hxed couplings Jij, the distribution of Q (with respect 
to thermal fluctuations) is related, at low temperature T, to the hrst gap 
between the two largest eigenvalues Ai^jv and A 2 ^n of the matrix J 


P{Q,N)r.-P 


^ P(^q) = e-i^(Ai,iv-A2.iv)(l--?'‘) _ 


( 12 ) 


Hence we see from (12) that the full distribution of the overlap in this model 


is directly related to the PDF of the hrst gap di^M = ~ A 2 ^n of Gaus¬ 

sian random matrices Jjj, the natural ensemble being here GOE (/3 = 1). 
Indeed, after averaging the distribution P{Q,N) in (12) over the random 
couplings Jjj, one obtains nothing else but the Laplace transform (with 
Laplace parameter oc (1 — q^)/T) of the PDF of the hrst gap, PGAPi^, N), 
studied here [see Eq. 


3. Two different scaling regimes: bulk and edge 

Gomputing the DOS /Odos(^) Af) for hnite N is, for a generic value of /3 in 
(Q, a very challenging task. This could be done for the special case /3 = 2 in 
Ref. m using the method of orthogonal polynomials, which unfortunately 
can not be extended to other values of /?. In spite of this difficulty, the 
main features of /Odos(^) Al), for large N, can be characterized by means of 
heuristic arguments. 

For this purpose, it is hrst useful to recall that the huctuations of the 
eigenvalues are characterized by two different scales depending on their loca¬ 
tion in the spectrum: (i) in the bulk for \i/y/N = 0(1) and |Aj| < y/2N and 
(ii) at the edge where [A* ± V‘^A[\ = 0{N~^A)^ The existence of these two 
scales manifests itself in various observables associated to the eigenvalues of 
the Gaussian /3-ensemble (Q, including their mean density dehned as: 
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One has obviously p(X,N) = p{—X,N) and one can further show that, for 
large N, it exhibits two distinct regimes (for A > 0) [T8l [l9l |20l [2T] 


p{X, N) ~ < 


1 


^/N 

V2 

liV5/6 


Pbulk 


A 

Vn 


, X < & A = 0{Vn) , 


(14) 


Pedge ((A - \/^)\/2iVV6) , |A - Vm\ = 0{N-^/^) . 


In Eq. (14), /?bulk(a^) is the Wigner semi-circle |T8l IT^: 

1 
vr 


Pbuik(a:) = Pw{x) = - \/2 - , 

TT 

independently of /3, while /?edge(a^) is given by [20l EU [22], 

' [Ai'(x)]^ — xAi^(x) -I- ^Ai(x) (l — dtAi(t)) ,,0 = 1 

[Ai'(x)]^ — xAi^(x) , 


(15) 


Pedge(2^) ~ 


/3 = 2 


(16) 


.-1/2 


^[Ai'(Kx)]^ — KxAi^(Kx) — ^Ai(Kx) dfAi(f)^ , 


/? = 4 


where k = 2^/^. Its asymptotic behaviors are given by 


Pedge(,^j y X 
TT 


X —OC 


(17) 


ln(Pedge(a;)) ~ ’ X->CX) 


Interestingly, one can check that these two regimes for p[X, N), the “bulk” 
one and the “edge” one in Eq. (141, perfectly match when A approaches the 
value \/2iV from below. Indeed, when A —?■ y/2N from below, p{X, N) can 
be replaced by the Wigner semi-circle (15), which gives: 


p{x,N) - N-^/^(Vm-x) ,x^V^ . 

TT \ / 


1/2 


(18) 


This behavior (18) coincides with the left tail of the scaling function pedge(a^) 
in Eq. (17). Indeed, when the deviation from \/2N is large, y/2N — A ~ 
0{Vn), we can substitute in the second line of Eq. ( |l4| the left tail asymp¬ 
totic behavior of Pedgeix) in 0. which gives 

p{X, N) ~ V2N-^/^^ (y2y/^{V^ - A))A ^ Vm~ , (19) 
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which after a trivial rearrangement coincides with Eq. (18). Note also that 
the right tail of /?edge(a^) in Eq- 0 matches, as it should, with the right 
tail of the TW distribution for /5-ensemble |19) . 

Similarly to the density 0. one expects that, for large N, pnosif^^) 
exhibits two different scaling regimes (see figure]^: (i) a bulk regime, where 
r oc \/]V and (ii) an edge regime where r = They can thus be 

summarized as follows 


PDOsir,N) ~ < 


1 _ f r 

Pbulk 


y/N 


VnJ ’ 


c < r < 2\/2N , 


( 20 ) 


,v/2iV-5/6pedge(rv/2iVi/6) , ^ = o(A^-V6) 


for some real c < 2\/2 and where Pbulk(a^) and pedge(i^) are two different 
scaling functions. 

Let us first investigate the bulk regime. For r oc y/N in (j^, one expects 
that /Odos(^) is insensitive to the fluctuations of Amax which are of order 
around the value y/2N. Therefore, in Eq. t he PDF of Amax 
can be simply replaced by a delta function (5(Amax ~ v2iV). It thus follows 
that P'Dosi'TjN) ~ p{y/2N — r,N) where p{X,N) is simply the density of 
eigenvalues in (|13|). Therefore, from (14) together with (15) one expects 


that /9DOs(^)-Ai) takes the scaling form given in the first line of Eq. (20) 
where /5buik(3;) is a shifted Wigner semi-circle: 


/Obuik(a;) = pw{V2-x) = -Jx{2y/2-x) , 

TT ^ 


( 21 ) 


independently of (5. 

On the other hand, in the edge regime when r ~ the DOS 

will be sensitive to the fluctuations of Amax and we expect that Ped^ei^) is 
a non trivial function, as it was shown to be the case for /5 = 2 m (see 


also Eq. (24) below). Although we can not compute explicitly this scaling 
function /5edge(f^) for other values of /5 7 ^ 2 we can extract its asymptotic 
behaviors for both small and large arguments, which we discuss in the next 
section. 

In Ref. m these results were obtained by exact analytical calculations 
for /5 = 2. Here we have performed numerical simulations for different values 
of /5 = 1,2 and 4 using the tridiagonal random matrices representation [ 8 ]. In 
Fig. we show a plot of poos (r', N) as a function of r and different values 
of N that corroborates the scaling forms in Eq. (20) for three different 


values of /5 = 1 (red), /5 = 2 (green) and /5 = 4 (blue). The central panel 
indicates that, in the bulk regime, Pdos(^)-A^)) correctly rescaled, converges 
to a shifted Wigner law, independently of /5. This is in full agreement with 













r 




Figure 2. The different scaling regimes for pdos as described in Eq. (20l. Left 
panel: Plot of pdos(^j-^) as a function of r for three different values of N and 
three different values: /3 = 1 (red), /3 = 2 (green) and /3 = 4 (blue). Central 
panel: (bulk regime) Plot of the same quantities (after rescaling) '/N pdos (^j N) 
as a function of r/^/N and for different value of jS with the same color code as the 
one used in the left panel. Right panel: (edge regime) Plot of the same quantities 
(after yet another rescaling) N )/as a function of 


the scaling form in the first line of Eq. (20). Finally, the right panel shows a 


plot of /Jdos(^) for small r, which is in a good agreement with the scaling 
form given in the second line of Eq. (20). It also indicates that the limiting 


function pedge(i^) depends explicitly on /3. 

We conclude this section by mentioning that the typical fluctuations 
of the first gap are naturally expected to scale as = ^i,N — ^ 2 ,N ~ 
as the fluctuations of Ai^jv = Amax and A 2 ^n around \/2N are 
also of order [HE]. Hence, for large N we expect that PGAp(r, N) 

takes the scaling form: 


PGAp(r,iV) = \/2iVfo%p (r\/2iVfo6^ 


( 22 ) 


where the factor \/2 in the argument has been chosen here for convenience. 
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according to the choice made for /5 dos(^)-^) in second line of Eq. (20) 


This scaling form (22) together with the scaling function ptyp{f) were ob¬ 
tained exactly for (3 = 2. It is natural to expect that this scaling form (22) 


holds for any value of /? > 0, with a /3-dependent scaling function Ptyp{r), and 
below we give some heuristic arguments to obtain the asymptotic behaviors 
of the scaling function ptyp (r) ■ 


4. Asymptotic behaviors of the scaling functions Pedge(i^) 

and Ptyp{r) 

We first begin to analyze the small f behavior of Pedge(f^)- From the 
definition of the DOS in Eq. (§ it is clear that its small f behavior is 
directly related to the probability that two eigenvalues - namely the first one 
and the second one - become extremely close to each other. From the joint 
PDF of the eigenvalues (Q, this probability vanishes as as a consequence 
of the short distance repulsion between eigenvalues, which comes from the 
Vandermonde term Wi j |Aj — Xj\^ in the joint PDF (j^. One thus expects 
that Pedgei^) ~ with an a priori unknown constant ap. 

On the other hand, it is reasonable to assume that there is a smooth 
matching between the edge region and the bulk region described by the 


shifted Wigner semi-circle law (21), as it is the case for the density of eigen¬ 


values [see the discussion between Eqs. (18) and (19)]. This means that the 
large f behavior of pedge(f^) has to coincide with the small argument of the 
shifted Wigner semi circle (21), which does not depend on (3. Hence one 
deduces that /5edge(f^) ~ for all (3. The asymptotic behaviors of the 

DOS in the edge regime can thus be summarized as follows: 


/5edge(f^) 


~ 


0/3 


+ o{f^) , f —>• 0 , 


(23) 


—^ -I- o(f^/^) , f —^ oo. 


TT 


lation of /Oedge(^) [TT] : 

21/3 


These behaviors (23) have been confirmed, for /3 = 2, by an exact calcu¬ 


Pedge(T) — 


TT 


f 


f{r,uf - 


q{u)f{f, u)du 


F2{x)dx , (24) 


where J- 2 {x) is the Tracy-Widom distribution associated to GUE [T] and 
q{x) is the Hastings-Mc Leod solution of the Painleve II equation 


J 2 {x) = exp 


poo 

\ 

— {u — x)q^{u)(lu 

J X 

’ 1 


q''{x) = 2q^{x) +xq{x) , 
q{x) Ai(x) for X —)> 00 , 


(25) 
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Figures. Plot of Pedgei^) for three different values of /3 (/3 = 1 in red, /3 = 2 in 
green and /3 = 4 in blue). Left panel: the data are plotted on a log-log scale and 
they exhibit a small f behavior compatible with our heuristic arguments in (231. 
Right panel: plot of Pedge{f) in terms of r with the asymptotic behavior for large 
f compatible with our heuristic arguments (231. 


while f{r,x) satisfies 

dlRf, x)-[x + 2q^{x)]f{f, x) = x) , 

f{r,x) ~ 2“^/®y'7rAi(x — f) . 


(26) 


It was further shown in m that /(f, x) can be expressed in terms of the Lax 
pair associated to the Painleve XXXIV equation. From this exact expression 


(241 one can derive the asymptotic behaviors announced in (23) and compute 


explicitly the amplitude 02 = 1/2 |11) . 

In Fig. I^we show a plot of Pedge(L) as a function of f for three different 
values of /3 = 1,2 and 4 exhibiting the small f behavior (in the left panel) 
and the large f behavior (in the right panel). In each case, /5edge(L) was 
computed by sampling 2.10^ independent random matrices of size 100 x 100. 
Both panels show a good agreement with our predictions in (23). We also 
notice on the right panel of Fig. |^that pedge(L) exhibits clear oscillations 
for /? = 4, while the curves for /3 = 1 and 2 are much smoother. 
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We now turn to the small f asymptotic behavior of the gap distribution 
Ptyp(f^), which can be argued to coincide with the small r behavior of the 
DOS, Pedgei^)- Indeed, for small f, the contribution to pedge(f^) comes only 
from the gap between the two largest eigenvalues because the others (the 
third, the fourth etc. eigenvalues) are too far, as a consequence of the short 
range repulsion coming from the Vandermonde term in Eq. Q. Hence one 
expects that ptyp{f ) ~ ap with the same amplitude as above in (23). 


On the other hand, the large f behavior of Ptyp(f^) can be obtained, to 
leading order, by assuming that for large separation f the two first eigen¬ 
values A 2 ^n and = Amax become statistically independent. Hence one 
expects: 


PGAp{r,N) = / Pr.[A2,Ar = A,Amax = A-hr]dA 

J —OO 

/ oo 

Pr.[A2,Ar = A]Pr.[Amax = A r]dA 

-OO 

Pr.[Amax = r] . 


r\j 

r^oo 


(27) 


r\j 

r^oo 


Therefore, we expect that right tail of Ptyp{f) coincides with the right tail of 
the TW distribution which goes for large argument as InT'^(x) ~ — 

To summarize, we obtain the asymptotic behaviors of Ptypir) as: 


Ptyp{f) ~ a/3 -h o(f^) r 

Inptyp(f^) ~ + o(f^/^) , r 

O 


(28) 


-|-oo 


while the computation of ap as well as the subleading corrections, both for 
small and large arguments, for any (3 remains challenging. 

In Ref. we obtained an exact expression of ptyp^f) for the special 
case (3 = 2 as 


Ptyp{r) = 


21/3 


TT 


f 


Pi-r,x) - 


q{u)f{—f, tt)drx 


F 2 {x) dx . (29) 


From this exact formula (29) we could not only check the above asymptotic 


behaviors (28) but also obtain the sub-leading terms as: 


Ptyp{r) = < 


-I- a4f^ -|- 0(7®) 


A exp 


— f3/2 _|— 

3 3 


21 

"32 


1 _ H^f-3/4 + 0(^i-3/2) 


1536 


( 30 ) 
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Figure 4. Comparison of ptyp('r) for the different /3-ensemble with 2.10^ realizations 
for a 100 X 100 matrix for each case. (/3 = 1 in red, /3 = 2 in green and ,3 = 4 in 
blue) Left panel: log-log plot of Ptyp{r) with the asymptotic behavior for small 
r. Right panel: log plot of ptyp in terms of with the asymptotic behavior 
for large f. We also plot the asymptotic behavior for large f obtained in the (3 = 2 
case Eq. (IM I. 


where the first and second lines correspond respectively to the small and 
large f behaviors. In (30), the amplitude A is given by ^ 


TT, 


where Ci^) is the derivative of the Riemann zeta function, while 04 can be 
expressed in terms of integrals involving q{x) with the result 04 ~ —0.393575.... 
It should be noticed that a different expression, somehow more complicated, 
of the PDF of the first gap Ptyp(f^) had been obtained previously in Ref. |23| . 
involving also Painleve transcendents. It remains an open question to show 
that these two expressions do coincide. 


5. The density of states and the gap for a fixed value of Amax 

Up to now, we have studied the DOS /Odos(^)-^) and the PDF of the 
gap pgap(?') N) averaged over the value of Amax- It is also interesting to look 
at these quantities for fixed value of Amax = V, the corresponding quantities 
being denoted as and pQAp{r\y, N). We naturally expect the 
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scaling forms, valid for all /3 > 0: 

PDOs(r-|y, A^) = iV-5/6^ /iedge(V2iV^/®r|^/2A^V6(y _ ^)) ^ (31) 

PGAp{r\y, N) = ArV6^ p,yp{V2Ny^r\V2N^/^{y - V^)). (32) 

From the results of Ref. im one can compute explicitly these scaling func¬ 
tions as: 


21/3 

Pedge{r\x) = - T2{x) 

TT 


R{x) 


f + 


R{x) 

q^{x) 


P{f, x) - x)g{f, x) 

q{x) 



(33) 


with R{x) = fP q‘^{u)du and q{x)g{f,x) = —f JP q{u)f{r,u)du . Simi¬ 
larly, for the gap one finds: 


21/3 

Ptypir\x) = - R2,s 

TT 



2^fs9s + 
Qs 



(34) 


where we haved now used the shortand notations for the shifted quantities 
(hence the subscript s) qs = q{x — f), Rg = R{x_^ f), R 2 ,s = J^ 2 {x — f), 


fs = f{—r,x — f), gs = g{—f,x — r). In Fig. (|y we show a plot of the 
constrained quantities: pedge(i^|0) in the left panel and Ptyp(i^|0) in the right 
panel when Amax = V2N ie x = 0 for three different values of /3 (/3 = 1 
in red, /3 = 2 in green and /? = 4 in blue) computed by sampling 2.10^ 
independent tridiagonal random matrices of size 100 x 100 [8] for each /3 and 
only kept the events when |Amax ~ <0.1. Both panels show a good 

agreement with our predictions in 34) for /3 = 2 (green solid line) |11| . 


6. Conclusion 

To conclude, we have studied the phenomenon of “crowding” of the eigen¬ 
values near the largest eigenvalue Amax for random matrices belonging to the 
Gaussian /3-ensembles (Q. In particular we focused on the DOS ([^, which 
is the density of eigenvalues seen from the largest eigenvalue, and the PDF 
of the first gap, these two quantity being related for any finite N through 
the relation (Q. Based on exact results obtained for Hermitian matrices, 
corresponding to /3 = 2, obtained in m, and using heuristic arguments 
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Figures. Comparison of pedge(?^|0) (left panel) and ptyp(?^|0) (right panel) for 
three different values of /3 (/? = 1 in red, /3 = 2 in green and /3 = 4 in blue). We 
have performed 2.10^ realizations of a 200 x 200 matrix. For each case, we keep 
only events when lAmax ~ < 0.1 (dots). We also plot in solid line the exact 

result for /3 = 2 (331 and ( [M| obtained in [TT] (using the numerical values of q{x) 
and T 2 {x) from Ref. m)- 


we obtained a general description of these two quantities in the large N 
limit. We also presented results of numerical simulations supporting our 
arguments. As we have seen in section]^ these quantities and 

PQfi^p(r, N) enter naturally into the computation of physical observables in 
the fully connected spherical spin-glass model and it will be particularly in¬ 
teresting to explore further the implications of our results to the relaxational 
dynamics of this model. Another application of the techniques developed in 
which provided a detailed analysis of an orthogonal polynomial system 
initially introduced in |12) . concerns the level curvature distribution at the 
soft edge of random Hermitian matrices which also involves the same orthog¬ 
onal polynomials |25) . Finally, it will be interesting to extend the present 
study to other ensembles of random matrices, like the Laguerre-Wishart en¬ 
semble. Hence we hope that these results will motivate further studies of 
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near extreme eigenvalues. 

Acknowledgments. We acknowledge support by ANR grant 2011- 
BS04-013-01 WALKMAT and in part by the Indo-French Centre for the 
Promotion of Advanced Research under Project 4604 — 3. GS acknowledges 
support from Labex-PALM (Project Randmat). We also thank Y. V. Fyo¬ 
dorov for useful discussions. 


References 


[1] C. A. Tracy, H. Widom, Level-spacing distributions and the Airy kernel, Com- 
mun. Math. Phys. 159, 151 (1994). 

[2] C. A. Tracy, H. Widom, On orthogonal and symplectic matrix ensembles, Com- 
mun. Math. Phys. 177, 727 (1996). 

[3] G. Ben Arons, A. Dembo, A. Guionnet, Aging of spherical spin glasses, Probab. 
Theory Relat. Fields 120, 1 (2001). 

[4] D. S. Dean, S. N. Majumdar, Large deviations of extreme eigenvalues of random 
matrices, Phys. Rev. Lett. 97, 160201 (2006). 

[5] D. S. Dean, S. N. Majumdar, Extreme value statistics of eigenvalues of Gaussian 
random matrices, Phys. Rev. E 77, 041108 (2008). 

[ 6 ] S. N. Majumdar, M. Vergassola, Large deviations of the maximum eigenvalue 
for Wishart and Gaussian random matrices, Phys. Rev. Lett. 102, 060601 
(2009). 

[7] S. N. Majumdar, G. Schehr, Top eigenvalue of a random matrix: large devia¬ 
tions and third order phase transition, J. Stat. Mech. P01012 (2014). 

[ 8 ] 1. Dumitriu, A. Edelman, Matrix models for beta ensembles, J. Math. Phys. 43, 
5830 (2002). 

[9] S. Sabhapandit, S. N. Majumdar, Density of near-extreme events, Phys. Rev. 
Lett. 98(14), 140201 (2007). 

[10] A. Perret, A. Gomtet, S. N. Majumdar, G. Schehr, Near-Extreme Statistics of 
Brownian Motion, Phys. Rev. Lett. Ill, 240601 (2013). 

[11] A. Perret, G. Schehr, Near-Extreme Eigenvalues and the First Gap of Hermi- 
tian Random Matrices, J. Stat. Phys. 156:843-876 (2014). 

[12] G. Nadal, S. N. Majumdar, A simple derivation of the Tracy-Widom distribu¬ 
tion of the maximal eigenvalue of a Gaussian unitary random matrix, J. Stat. 
Mech., P04001 (2011). 

[13] Y. V. Fyodorov, private communication after a talk given by one of us at the 
conference Random Matrix Theory : Foundations and Applications, Krakow, 
(2014). 

[14] G. Monthus, T. Garel, Typical versus averaged overlap distribution in Spin- 
Glasses: Evidence for the droplet scaling theory, Phys. Rev. B 88 , 134204 (2013). 



16 


[15] L. F. Cugliandolo, D. S. Dean, Full dynamical solution for a spherical spin- 
glass model, J. Phys. A: Math. Gen. 28, 4213 (1995). 

[16] J. Kurchan, L. Laloux, Phase space geometry and slow dynamics, J. Phys. A: 
Math. Gen. 29, 1929 (1996). 

[17] Y. V. Fyodorov, A. Perret, G. Schehr, in preparation. 

[18] M. L. Mehta, Random Matrices, 2nd Edition, Academic Press (1991). 

[19] P. J. Forrester, Log-gases and random matrices, Princeton University Press, 
Princeton, NJ, (2010). 

[20] M. Bowick, E. Brezin, Universal scaling of the tail of the density of eigenvalues 
in random matrix models, Phys. Lett. B 268, 21 (1991). 

[21] P. J. Forrester, The spectrum edqe of random matrix ensembles, Nucl. Phys. 
B 402(3), 709 (1993). 

[22] P. J. Forrester, N. E. Frankel, T. M. Garoni, Asymptotic form of the density 
profile for Gaussian and Laguerre random matrix ensembles with orthogonal and 
symplectic symmetry, J. Math. Phys. 47(2), 023301 (2006). 

[23] N. S. Witte, F. Bornemann, P. J. Forrester, Joint distribution of the first and 
second eigenvalues at the soft edge of unitary ensembles. Nonlinearity 26, 1799 
(2013). 

[24] This was done using the tabulation of q{x), R{x) and J^ 2 ix) which can be 
found on the webpage of M. Prahofer: http://www-m5.ma.tum.de/KPZ 

[25] Y. V. Fyodorov, Level curvature distribution: from bulk to the soft edge of 
random Hermitian matrices, Acta Phys. Pol. A 120(6), 100 (2012). 



